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Optical trapping is an indispensable tool in physics and the life sciences. However, there is 
a clear trade off between the size of a particle to be trapped, its spatial confinement, and the 
intensities required. This is due to the decrease in optical response of smaller particles and the 
diffraction limit that governs the spatial variation of optical fields. It is thus highly desirable to 
find techniques that surpass these bounds. Recently, a number of experiments using nanophotonic 
cavities have observed a qualitatively different trapping mechanism described as “self-induced back- 
action trapping” (SIBA). In these systems, the particle motion couples to the resonance frequency of 
the cavity, which results in a strong interplay between the intra-cavity field intensity and the forces 
exerted. Here, we provide a theoretical description that for the first time captures the remarkable 
range of consequences. In particular, we show that SIBA can be exploited to yield dynamic reshaping 
of trap potentials, strongly sub-wavelength trap features, and significant reduction of intensities seen 
by the particle, which should have important implications for future trapping technologies. 


Optical trapping is one of the most important ex¬ 
perimental tools in physics and life sciences because it 
enables precise control over small dielectric particles [1]. 
Famous examples of its use are optical levitation and 
cooling of nanoscale particles [2-6], trapping of bacteria 
[7] and cells [8], optical sorting in microfluidic channels 
[9], the manipulation and stretching of DNA [10], and 
recently, even trapping of individual HIV-1 viruses 
[11]. However, the difficulty of trapping a particle 
generally increases with decreasing size, due to the 
decreased optical response of the particle. This requires 
a commensurate increase in field intensity to maintain 
trap stability, and leads to associated problems such as 
thermal or material damage. Another limiting factor is 
the diffraction limit, which constrains the length scale 
over which helds can vary, and thus the stiffness or 
possible spatial features that a trap can possess. 

A number of experiments in recent years have mi¬ 
grated from trapping in free-space beams to the helds 
generated in nano-optical resonators [12-17] as illus¬ 
trated in Fig. 1. Such a paradigm can enable some 
technical advantages. For example, the resonator 
allows one to build up a higher intensity seen by the 
particle within the structure compared to the input, 
thus relaxing input power requirements. Engineering 
the nanophotonic structure also provides some hexibility 
over the held prohle, and thus the trapping potential. 
However, it is clear that simply replacing the input held 
with the enhanced one does not relax any requirements 
from the standpoint of intensity seen by the particle. 
Therefore it remains an open question whether one 
can circumvent these seemingly fundamental trade 
offs between particle size and the intensities required 
to achieve given trap depths, frequencies, and spatial 
conhnement. At the same time, doing so would have 
signihcant implications for optical manipulation as a 


tool in physics, chemistry and biology. 



FIG. 1: Schematic illustration of trap configurations. 

a) A dielectric particle trapped with an optical tweezer in free 
space, b) a plasmonic cavity {e.g., a metallic bowtie antenna), 
c) a photonic crystal cavity. 


In this context, a number of experiments have observed 
qualitatively new trapping behavior in nanophotonic 
cavities [12, 13]. The key physics is that the position of 
the trapped particle alters the resonance frequency. This 
results in a “self-induced back-action” (SIBA) effect in 
which the motion dynamically affects the build up of 
intra-cavity intensity, and thus the optical force exerted. 
However, the involved trapping mechanism and its range 
of consequences has hardly been explored. 

Here, we develop a general theoretical model for 
SIBA. Using such a model, we show how parameters 
can be chosen to maximize the effects of back-action, 
and that a single “back-action parameter” rj cx 
proportional to the resonator quality factor and the ratio 
of particle to cavity mode volumes, characterizes the 
performance of any optimized system. In particular, the 
back-action parameter indicates how many line widths 
the particle can shift the cavity resonance frequency due 
to its movement. For large rj, large shifts in the cavity 
detuning relative to the laser frequency as the particle 
moves can induce strong changes in the intra-cavity 
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intensity. Under these circumstances, and when properly 
optimized, such a trap yields very different trade-offs 
between intensities, trap depth, and confinement, which 
should have significant consequences for optical trapping 
technology. Specifically, we show that back-action can be 
exploited to create traps with strongly sub-wavelength 
spatial features, even if the cavity mode itself obeys 
the diffraction limit. The spatial features of the trap 
can also be dynamically shaped using only changes in 
laser frequency. Furthermore, the particle can effectively 
be trapped in a dynamical intensity minimum, even 
if it is nominally high-intensity seeking, which can 
strongly reduce the effects of photo-thermal damage. 
Finally, we discuss the possibilities for implementation in 
nano-plasmonic (Fig. lb) and photonic crystal (Fig. Ic) 
systems. 

MODEL OF TRAPPING IN NANOSCALE 
RESONATORS 

We first briefly review the properties and limits of trap¬ 
ping with free-space optical tweezers. Subsequently we 
will present our model of trapping in nano-resonators. 
Considering a small dielectric particle whose dimensions 
are much smaller than the optical wavelength c? ^ A its 
response is that of a point dipole with induced dipole mo¬ 
ment pind = a(uj)E{x). The well-known time averaged 
potential for optical trapping in free-space, such as by an 
optical tweezer (see Fig. la) reads [1]: 

UT{x) = -\Re{a{u:))\E{x)\^ ( 1 ) 

where a(w) is the frequency dependent polarizability 
of the particle and \E{x)\^ oc I{x) is the peak electric 
field amplitude squared at the particle position x and 
at frequency w, which is proportional to the intensity 
I{x). In the following, we will focus on the case where 
the polarizability is positive and largely frequency inde¬ 
pendent, which models well a typical dielectric particle. 
In this case, the dielectric particle is trapped around 
points of local maximum intensity. For sub-wavelength 
particles, the polarizability is proportional to particle 
volume, a(uj) oc V. It can thus be seen that the trapping 
of smaller particles requires a commensurate increase in 
intensity to maintain a fixed trap depth. Furthermore, 
the spring constant around the trap minimum Xmin, 
^spring = UT(a;min) ^ , in addition to being pro¬ 

portional to the beam intensity / and particle volume 
V, is at best proportional to the inverse of the optical 
wavelength squared, as the diffraction limit sets this as 
the minimum scale over which free-space optical fields 
can vary. 

We now examine the case where the particle is 
trapped in a nanoscale cavity. Our formalism is quite 


general, covering equally systems such as plasmonic and 
photonic crystal cavities, and trapping in vacuum or 
fluid environments. Qualitatively, the new feature of 
such a system is that the resonance frequency of the 
cavity depends on the particle position, enabling the 
particle motion to feed back on its trapping potential. 
We then distinguish the regimes in which this system 
gives rise to standard optical trapping as in Eq. (I), 
versus a novel “back-action” trapping mechanism. 

A general model of this system is given by follow¬ 
ing Hamiltonian: 

H = Hu}Jx„)a''' a + HJE^Eq -I- H - 

' ' 2m 

( 2 ) 

where ujc{xp) is the resonance frequency of the optical 

cavity as a function of particle position Xp and a is 

the annihilation operator of the cavity mode. Kex 

denotes the decay rate of the cavity into some particular 

external channel (such as free-space radiation, coupling 

fiber, etc.), which also serves as the source of injection 

of photons into the cavity with number flux Eq and 

2 

frequency ojl- The last term, Akin = describes the 
kinetic energy of a particle with momentum p and mass 
m. In addition to the external coupling, we assume 
that the cavity has an intrinsic loss rate Kin, such as 
through material absorption or scattering losses. The 
total cavity linewidth is thus k = Kin + Kex- In principle, 
the particle also contributes a position-dependent loss 
term K{xp) due to its scattering of light out of the cavity 
mode. While this term could be explicitly included in 
the analysis, this position-dependent effect is negligible 
under reasonable conditions as the scattering rate oc ^ 
rapidly falls off for sub-wavelength sizes, as shown in 
the Appendix. Thus, the quality factor of the resonator 
is defined as Q = where uic is the empty cavity 
resonance frequency. 

The system dynamics under the Hamiltonian of 
Eq. (2) and system losses are described by standard 
Heisenberg-Langevin equations [18]. As the regime of 
interest for trapping is far from any quantum behavior, 
we proceed to solve their classical expectation values. 
We neglect damping of the mechanical motion and the 
effect of a thermal environment, which do not influence 
the optical force and can be added independently later 
on. The equations of motion then read 

dx^ _ p_ . , 

dt m 

^ = -n{xp)hw'^{xp) ( 4 ) 

= i(a;L - Wc(a;p))^ - |/3-b Vk^Eo (5) 

where j3 = {a) is the expectation value of the photon 
amplitude while n = |/3p is the expectation value 
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of the photon number in the resonator. We note 
that even in state-of-the-art photonic crystal cavities, 
the achievable quality factor Q = ^ « 10® results 
in decay times of ~ 1 ns that are significantly 
shorter than the timescales of motion [19]. Thus, this 
motivates an approximation ^ ~ 0 where the cavity 
is able to instantaneously respond to the particle motion. 

Before solving equations (3)-(5), we want to exam¬ 
ine how strongly the particle affects the resonance 
frequency. To quantify this, we compare the frequency 
shift ScOcixp) = ujc{xp) — uJc with half of the line width 
^ of the resonator, where uJc is the resonance frequency 
of the empty cavity. Within lowest order perturbation 
theory, where the particle induces a frequency shift 
much smaller than the bare cavity frequency, it can be 
shown that (Appendix) 

2Sujc(xp) . , , . 

- — = -vf{x). (6) 

K 

Here, we have defined the dimensionless back-action 
parameter rj = Q, where Vm is the cavity mode 
volume. f{x) is the dimensionless spatial intensity 
profile of the empty cavity, normalized to be 1 at 
the intensity maximum. Thus, as 0 < f{x) < 1, 
the back-action parameter ry characterizes how many 
linewidths (half-width half-maxima) the particle can 
shift the resonance frequency of the cavity moving 
from the minimum f{x) = 0 to the maximum of the 
mode profile. For sub-wavelength dielectric particles the 
polarizability a{u!) oc cqV is proportional to the particle 
volume, with the pre-factor depending on the particle 
refractive index and shape [20]. Thus, achieving a large 
back-action parameter requires a sufficient combination 
of large cavity quality factor and ratio of particle to 
cavity mode volume, ry oc ^ • Q. When the particle size 
is larger than kr > Q~^ (with k = 27r/A), the effect on 
the quality factor due to light scattering by the particle 
cannot be neglected anymore, and this case is considered 
further in the SI. 


The expectation value of the intra-cavity photon 
number n{xp) = J/!]^ reads: 


( \ 4ii/Q Kex 


i-H 


(vfixp) 


+ A 


(7) 


where we have defined the dimensionless detuning be¬ 
tween the laser and empty cavity frequencies, A = 
^ From Eq. (7), one sees that there are certain 
positions of the particle Xr that cause the driving laser 
to become resonant with the (frequency-shifted) cavity, 
A + r]f(xr) = 0, and where the intra-cavity photon num¬ 
ber is maximized. We call these positions the resonant 
positions Xr, which can be chosen by adjusting the laser 


frequency Note that in arbitrary dimensions the res¬ 
onant positions Xr are contour points/lines/surfaces in 
1D/2D/3D and follow the symmetry of the mode profile, 
see Fig. 2. Inserting Eq. (7) into Eq. (4) and integrating 



FIG. 2: Back-action trapping in the fundamental mode 
of a Fabry-Perot cavity, with dimensionless intensity pro¬ 
file f{x) — cos^(fcx) (dashed blue curve), a) In the regime 
of small back-action parameter (?y = 0.1), the intra-cavity 
intensity is not significantly affected by the particle motion. 
Thus, the local intensity I{Xp) seen by the particle (red) is 
directly proportional to f{x), while the trapping potential 
U{x) oc —f{x) (yellow). For an increasing back-action param¬ 
eter ry 3> 1 the seen local intensity I{xp) forms sharp peaks 
centered around the resonant points Xr and the trapping po¬ 
tential U{x) converges to a square well potential, b) Spectra 
of intra-cavity photon number n{ujc{xp)) taken at the instan¬ 
taneous particle positions shown in the ?y = 10 and ?y = 50 
cases, respectively. For jy = 10, we consider the case where 
the particle is instantaneously located at one of the resonant 
positions x,.i, such that the laser frequency is resonant with 
the cavity at this moment to generate a large intra-cavity 
intensity. For ly = 50, the particle is far from the resonant 
positions, and the large detuning of the laser from resonance 
strongly suppresses intra-cavity intensity. The vertical scales 
of these plots are in arbitrary units. 


the negative force with respect to Xp yields the general 
potential for trapping in resonators: 


U{x) = —2hEQ^^ arctan 


Vfix) + A 


( 8 ) 


We will proceed by looking at different regimes of this po¬ 
tential: First we consider the regime where the particle 
induces a shift on the cavity resonance frequency that is 
negligible compared to its linewidth, which corresponds 
to 7y ^ 1 from our definition in Eq. (6). Then, the move¬ 
ment of the particle does not significantly change the 
intra-cavity intensity, which recovers the optical tweezer 
regime. In particular, expanding Eq. (8) for small rj, 
one finds that Ut{x) = —2H^Eq ^ fix). Using 

the definition of ry = and identifying \E{x)\^ = 

S-A^^QEq ^^^ 2 f{x) as the time averaged intra-cavity 
field amplitude, we see that Ut{x) reduces to the optical 
dipole potential in Eq. (I). In this regime, the potential 
depth increases linearly with Q (i.e., with rj), reflecting 
the effect of a built-up intra-cavity intensity. The differ¬ 
ent regimes are illustrated in Fig. 2 where we choose the 
first harmonic of a Fabry-Perot cavity as a mode profile. 
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TRAPPING WITH BACK-ACTION 


We now investigate the very different trap properties 
that emerge in the regime 77 ^ 1 . 

An increase in the quality factor initially produces 
an increased trap depth for values Q < (at which 

point ry ~ 1). For larger values, however, rj ^ 1 and 
the arctan in Eq. 8 saturates between the values of ±|, 
yielding a trap depth of 6U = 27r/iEg^. Significantly 
for ?7 1, the depth no longer depends on Q nor the 

particle properties, and is only dependent upon the 
input intensity. The origin of this saturation can be 
understood by first considering Fig. 2, which shows 
that the intra-cavity intensity as a function of particle 
position forms sharp peaks around the resonant positions 
Xr for r] ^ 1. From Eq. (7) it follows that their width is 
in good approximation « ^ and it is only within 

this narrow spatial region (scaling like 77 “^ oc < 5 “^) that 
the cavity exerts significant forces on the particle. At 
the same time, the peak intra-cavity photon number 
at Xp = Xr (and thus the peak force) grows linearly 
with Q. Thus, the maximum work that the cavity can 
do to keep the particle in the trap, as a product of 
force and distance, becomes independent of Q in the 
high-back-action limit. 

Note that the trapping potential turns into an ap¬ 
proximate square well if the distance between the 
intra-cavity intensity peaks is larger than their width 
d = \xt :2 — Xri\ rif\x ) ' wells are (symmetrically) 
centered around the mode profile maximum Xq, see 
Fig. 2. Remarkably, the resonant positions Xr can be 
changed with laser frequency, which provides a conve¬ 
nient mechanism for dynamic trap shaping in contrast 
with conventional optical tweezers. 

Another interesting property of the trap in the high 
back-action regime is that around the minimum Xq of 
the potential, the intra-cavity photon number is strongly 
suppressed due to the large detuning from resonance. 
Thus, the particle is effectively trapped in a dynamical 
intensity minimum, despite the fact that it has positive 
polarizability and is thus nominally high-intensity seek¬ 
ing. This would have tremendous consequences in the 
reduction of thermal damage due to optical absorption 
by the particle. Motivated by this observation, we seek 
to quantify how much the time-averaged intensity seen 
by the particle can be reduced. 

We define the time-averaged experienced intensity 
{Iexp)t as the local intensity experienced by the particle 
at its position, averaged over one motional period T. It 


is thus given by 

{Iexp)t = ^ n{xp{t))f{xp{t))dt (9) 

where Xp{t) is a solution to the differential Eq. (4) 
together with Eq. (7). In order to proceed further, we 
consider a simple case of the fundamental mode of a 
ID Fabry-Perot cavity, f{x) = cos^{kx) with k = ^, 
where L = A/2 is the cavity length. Although we have 
switched to a specific model to illustrate the back-action 
mechanism, we believe the overall conclusions are 
generally valid. A finite temperature of the environment 
can be taken into account by averaging the results for 
different maximal kinetic energies Akin (kinetic energy 
of the particle in the trap minimum) according to a 
Boltzmann distribution. 


We have evaluated Eq. (9) by numerically solving 
the equations of motion (3)-(5). In Fig. 3, we plot the 
time-averaged experienced intensity {Iexp{v))t normal¬ 
ized by the value in the optical tweezer regime (/exp.r)*, 
as a function of back-action parameter 77 . As seen before, 
the optical tweezer regime is reached by taking 77 ^ 1 . 
To make a fair comparison, we enforce that the trap 
depths in the two cases are equal, SU{ri) = 6Ut- For 
a fixed Xr, the figure shows a significant reduction in 
time-averaged intensity for high back-action parameter, 
which also depends on the ratio of kinetic energy Ekin 
to trap depth 6U. In the high back-action regime, it is 
possible to derive an analytic expression (Appendix): 


lim (/exp(?7))t 

rf-^-oo 


2ceo f{Xr^ Ekin 
a(uj) \rixr)\ Xr 


( 10 ) 


A new feature of the back-action trap is the gradual 
decoupling between trap depth and the spatial region 
6x = \xt 2 — Xti\ (xti and Xt 2 are the classical turning 
points) to which the particle is confined. For large 
enough 77 they decouple completely since the classical 
turning points converge to the resonant positions (i.e., 
the edges of the square well) and thus 6x —> d. In this 
regime, confinement only depends on laser frequency, 
whereas trap depth only depends on laser power. This 
independent control again highlights the ability to 
dynamically reshape the trap. In contrast, in the optical 
tweezer regime, the trap depth, kinetic energy and 
confinement are inevitably connected. 


Instead of comparing the experienced intensity at 
fixed trap depth, we can also investigate the trade-off 
between intensity and confinement Sx = d in the large 
back-action limit. The locations of the trapping wells 
are always centered around the mode profile maximum 
xo = 0. For small Xr, an asymptotic expansion yields 
~ 2 khr- Thus, for high back-action and strong 
confinement, we obtain {Iexp)t « IfSy ^kin- In¬ 
terestingly, expanding Eq. (1) for the optical tweezer 
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around the bottom of a standing wave potential also 
produces (/exp,T)t « I{xo) « IfSy -^kin, which 
seems to indicate that no improvement is gained in 
intensity vs. confinement with back-action. 

Looking at Eq. (10), in the strong back-action regime, 
one of the factors of originates simply from the 
time T (X Sx that the particle takes to travel between 
the walls of the square well. This part of the scaling 
seems fundamental and cannot be improved within this 
model. On the other hand, the second factor of 
clearly originates from the vanishing of back-action 
effects around the maximum of the mode profile, as the 
frequency shift becomes insensitive to first-order changes 
in the particle displacement, f'{xo) = 0. We show that 
this factor is not fundamental, and can be eliminated by 
properly driving a second optical mode of the system. 


(fexp {v))t 

(lexp.T ) t 



FIG. 3: Time-averaged experienced intensity of a 
trapped particle. We plot the time-averaged experienced 
intensity as a function of back-action parameter (/exp(?7))t, 
normalized with the value in the optical tweezer regime rj 1. 
The two cases are set to have equal trap depth. The plot is 
numerically calculated for the case of trapping in the funda¬ 
mental mode of a Fabry-Perot cavity f{x) = cos^{kx) with 
resonant positions kxr = The back-action regime can en¬ 
able much lower average local intensities than in the optical 
tweezer regime. 


TWO MODE BACK-ACTION 

In this section we show how the scaling between 
experienced intensity and confinement can be improved 
to (/exp)t oc by using two different cavity modes for 
trapping. In order to obtain concrete results, we consider 
the simple geometry where the two modes consist of the 
first and second harmonics of a Fabry-Perot (see Fig. 4), 
although we believe that the conclusions hold quite gen¬ 
erally. We assume that each mode can be driven with its 
own laser, with amplitude Eqi and frequency wl^. As the 


equation for the intra-cavity fields /3i (generalized from 
Eq. (5)) of each mode are decoupled from one another, 
they can be separately integrated as in the single-mode 
case. Thus, the total potential Utot{x) = J2i=i 
is the incoherent sum of the potentials in Eq. (9) for 
each mode. To understand the relevant physics, it is 
sufficient to assume that the mode driving amplitudes 
i?oi, decay rates Kex, and back-action parameters are 
identical, although the concepts can be easily generalized. 

The interesting regime will be when the resonant 
positions of each mode are tuned by their respective 
driving laser frequencies such that each mode is respon¬ 
sible for providing one trapping wall. This is illustrated 
in Fig. 4b, where the left and right walls Xri and 
Xr 2 originate from the first and second cavity modes, 
respectively. Significantly, the well can be located far 
from the nodes/antinodes /((x) = 0 where the effects 
of back-action would vanish for either mode. In the 
following we will distinguish three different regimes 
concerning the ratio between the distance d = jxri — x-c 2 \ 
and the width ~ ^ of these intensity peaks illustrated 
in Fig. 4. 

We start by examining the high back-action regime, 
when the distance of the intensity peaks is much larger 
than their width, kd = k\xri — Xr 2 \ such that we 

encounter an almost perfect square-well potential as 
shown in Fig. 4b. It is straightforward to generalize the 
high back-action limit of Eq. (10) in the single mode 
case. As the particle is trapped far from points where 
back-action effects vanish (/((a;) = 0), we recover the 
improved scaling between experienced intensity and 
confinement, (/exp)t cx as already anticipated. 

In Fig. 5, we have illustrated the results of experi¬ 
enced intensity vs. confinement from full numerical 
simulations of Eqs. (3)-(5) (generalized to two modes). 
Here, the different points for a fixed back-action param¬ 
eter r] are obtained by variation of the input powers and 
resonant positions Xr (via the laser frequencies). Tuning 
the resonant positions to reduce d = \x-ci — Xr 2 \ indeed 
enables one to saturate the scaling of (/exp)t cx ^ as 
long as kSx > as illustrated in Fig. 5b. 

For kSx the optimal scaling seen in the nu¬ 
merics goes like (lexp.hb)* «: ■ The scaling with 

5x~^ resembles the optical tweezer case, but the intensity 
is suppressed by a factor of rj. We call this the “harmonic 
back-action regime” (see Fig. 4a). To understand this 
case, we first note that the particle moves by a small 
enough amount around the trap minimum that the 
forces from each mode can be linearized around small 
displacements to yield a harmonic trap. Furthermore, for 
small displacements, the total time averaged experienced 
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intensity (/exp)t = Ej(-^exp,*)t ~ Y.i hi^o) is just the 
sum of the intensities of the respective mode at the trap 
minimum Xp = xg. The associated spring constant is: 

fcopt = -J^'(xo) = X] + ni(xo)uj'Ji(xo) 

i 

( 11 ) 

where the sum goes over all trapping modes. The first 



breakdown of the dipole approximation in which all of 
these expressions are derived. It should be noted that 
an analogous “optical spring” effect has been reported in 
optomechanical systems [21, 22], where an optical cavity 
can exert large restoring forces for small displacements 
of a mechanical system. Here, the stiffness of the me¬ 
chanical mode itself plays the role of our second optical 
mode, and serves to keep the equilibrium position at a 
point of non-vanishing back-action {f'{xo) ^ 0) [21]. 

Exploiting the notion of a reduced wavelength, in 
the harmonic back-action regime one can immediately 
conclude that the scaling for average experienced in¬ 
tensity improves from (/exp,T)t oc {kSxy optical 

tweezer to (/exp.hb)i oc A more detailed opti- 

mization of the system shown in Fig. 3 and explained in 
the SI reveals that: 


('fexp,hb)£ 4 

(4exp,T)£ V 


(13) 


FIG. 4: SIBA with two optical modes, illustrated here for 
the first two modes of a Fabry-Perot cavity. Top: the mode 
profiles are given by fi{x) = cos^{kx), f 2 {x) = sm^{2kx). 
Green: the resulting optical trapping potential [/(a;). Bottom: 
intra-cavity intensities I{x) as a function of particle position. 
a) In the harmonic back-action regime, the distance between 
the resonant points is comparable to the width of the intensity 
peaks, fed ~ A 6) In the high back-action regime, the distance 
significantly exceeds the width, fcd 3> 

term n[[xg) is a new contribution to the optical spring 
constant fcopt originating from the change in photon num¬ 
ber with particle position around the trap minimum. In¬ 
tuitively, this back-action contribution to the spring con¬ 
stant is maximized by ensuring the photon number of 
each mode maximally changes around xg. This is roughly 
optimized by setting kd ^ such that xg corresponds 
to sitting half a cavity linewidth away from the resonant 
position Xr- Such an optimization yields (Appendix): 

kopt,z = ^^(Jexp,j)t „ 1 - fiixg)] . 

ceg Ji[xg) 

( 12 ) 

The first term in the brackets originates from the 
change in photon number with particle position, whereas 
the second term reduces to the optical tweezer spring 
constant given by Eq. (1): kx = U"{xg). Since 
f''{xg) ^ f'{xo)^ ^ it can be seen that the back- 
action contribution is a factor of rj larger. We can 
equivalently interpret this contribution as arising from 
an effective reduced wavelength Aeff , which enables 
the generation of trap features far below the diffraction 
limit. We emphasize that this effect originates from 
the rapid change in intra-cavity photon number with 
particle displacement rather than a change in the 
spatial mode itself (see Eq. (11)), and thus there is no 


for equal confinement and kinetic energy. We want to 
emphasize that to reach this optimal scaling, one should 
fix kd = k\x^i — Xr 2 \ ^ p- l'^ llil® way, one stays on 
the dashed line scaling shown in Fig. 5 and one achieves 
smaller confinement by turning up laser intensity while 
still maintaining the full back-action advantage. In con¬ 
trast, Fig. 5 also shows that by decreasing the distance 
between the resonant positions, kd the scaling devi¬ 
ates back towards the optical tweezer limit and the ben¬ 
efits of back-action vanish. 


CONCLUSION 

There have already been two types of systems, plas- 
monic cavities [23] and photonic crystal cavities [19], 
where SIBA has already been observed, and we now dis¬ 
cuss the potential s of merit associated with each. As 
the plasmon resonances associated with small metallic 
systems do not obey a diffraction limit, they are able 
to achieve strongly sub-wavelength mode volumes. On 
the other hand, realistic quality factors are limited to 
Q < 10 — 10^. At the same time, an upper bound on the 
validity of our calculation is that the particle size V < Vm 
does not exceed the mode volume, and thus we anticipate 
maximum possible values of ry ~ 10 — 10^ for such sys¬ 
tems. In photonic crystal cavities, the mode volume is 
limited by the diffraction limit to Vm ^ (§)^, while ex¬ 
tremely high quality factors of Q ~ 10® are possible [19]. 
This yields rj ~ 10,100,400 for a dielectric sphere with 
radius r ^ 6.5nm,15nm, 28nm (Appendix). There has 
been significant activity in recent years to develop design 
principles in order to tailor the spatial modes of plas- 
monic [23] and photonic crystal structures [24] for trap¬ 
ping. Combined with the potentially large back-action 
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FIG. 5: Time-averaged experienced intensity vs. confinement for two optical modes, a) Time-averaged experienced 
intensity in units of as a function of confinement kSx = k\xt 2 — Xti\- The individual points originate from different 

combinations of back-action parameter, laser power and detunings. The solid lines indicate the scalings in the optical tweezer 
regime and high back-action regime, where kd — k\xr 2 — ®ri| ^ The dashed line shows the optimized harmonic back-action 
regime, where kd ^ b) Illustration of figure 5a for a fixed rj and ifkin, and a schematic of the protocol to saturate the 
scaling bounds. The green arrows denote a decrease in distance between the resonant positions d = \xii — Xr 2 \, while the red 
arrows denote an increase in laser power. The ratio of filled blue area to whole area of the parabolae/square wells indicates 
the ratio of p^^). 1) high back-action regime: d and SU (laser power) decouple, the potential forms a square well 

and Sx —> d. 2) a significant increase in laser power prevents the particle from exploring the full square well potential. 3) 
decreasing kd < ^jr] suppresses back-action as the particle motion no longer shifts the cavity mode frequencies. 4) beginning 
of the optimized harmonic back-action regime. Maximum confinement at this point is achieved when « |. Maintaining 
the harmonic back-action scaling (dashed line) is achieved by increasing laser power. 


parameters achievable, we anticipate that our work will 
open up significant new opportunities for optical trap¬ 
ping. Finally, it would also be interesting to explore the 
use of large back-action parameter in other functionali¬ 
ties, such as particle detection and feedback cooling. 
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Here we provide additional supporting calculations that were not included in the main text. 


1 FREQUENCY SHIFT 

For small frequency shifts Sujdxp) = ajc(xp) — ojc com¬ 
pared to the unaltered resonance frequency Wc of the cav¬ 
ity, we can obtain Sujc{xp) from electromagnetic pertur¬ 
bation theory [1]: 


with Vm being the mode volume, c the photon veloc¬ 
ity, and (Tscat = the Rayleigh scattering 

cross-section (where /c = ^ is the wavevector of the 
incident light). We begin by comparing the relative 
effects of the position dependent scattering rate and 
cavity frequency shift on the intra-cavity photon number. 


ScJci^Xp'j 


Vc f d?rP(r) ■ E{r) 
^/dV^(e(a;,r)cu)|E(r)|^ 


(SI) 


where E{r) is the electric field of the empty resonator, 
e(w, r) is the dielectric function of the empty resonator, 
and P{r) is the additional polarization due to the pres¬ 
ence of the particle. If we take the particle to be small 
compared to the wavelength of the laser, the electric field 
across the particle is approximately constant and its re¬ 
sponse is equivalent to a point dipole with polarizability 
a{uj). As an example, for a dielectric sphere of volume V 
and refractive index n in vacuum, the polarizability can 
be exactly calculated, a(w) = For a given 

polarizability, one finds 


SoJc{Xp) 


tOcCcjuj) 

2eoVm 


fix) 


(S2) 


where f{x) = 




r. f{x) describes the 


maxx (e{oj,x)io)\E{x)\'^ 

dimensionless spatial intensity profile of the empty cavity, 
normalized to be 1 at the intensity maximum. Vm is the 
mode volume of the empty resonator and is defined as 
follows: 


I {e{uj,x)uj) \E{x)\^ 
maxa; ^ (e(w,x)w) \E{x)f 


2 SCATTERING RATE OF THE TRAPPED 
PARTICLE 


Fig. 2b in the main text shows how the particle 
motion shifts the resonance peak of the intra-cavity pho¬ 
ton number spectrum. In contrast, a position dependent 
scattering rate does not shift the peak, but instead alters 
its width and height. With this picture in mind we can 
neglect the effect of the position dependent scattering 
rate, if the change in scattering rate 5k(Xp) = Kscatixp) 
is much smaller than the frequency shift StOcixp) induced 
by the same particle movement. Using Eq. (S4) and 
Eq. (S2) and comparing these two quantities yields: 


l^^scat {Xp) I 

\SuJc{Xp)\ 


(kr)^ <C I, 


(S5) 


which allows us to neglect the position dependence of 
the scattering rate for sub-wavelength particles. 


Nonetheless we have to consider the reduced qual¬ 
ity factor of the resonator-particle system due to 
scattering of light. The total cavity decay rate is 
K = ^ + Kscat(a;o) where Q = is the quality 

factor of the empty cavity. Thus, the back-action 
parameter reduces to 


V = Q 


a(uj) 
eoVm 1 -p 


I 

Kex-I-Kint 


(S6) 


From Eq. (S4), assuming that a{u!) « eoV and writing 
Vm = V , where v tells us how close the light is fo¬ 
cused to the diffraction limit, the scattering rate reads 


Here we consider the scattering of light by the trapped 
object itself, which decreases the cavity quality factor by 
contributing to its loss rate n = Kex + Kint + Kscat(a^p)- 
The scattering rate for sub-wavelength particles reads: 


^scat(^p) — ^scatC 


fjXp) 

Vm ’ 


(S4) 


^scat(Xp) ~ ^scat(^o) ~ 9 (kv) OJl,. i^P 

/ ITT E 

Inserting this into Eq. (S6) finally yields 


4 Q (kr)^ 


(S8) 



2 


and is plotted in Fig. 1 and Fig. 2 for = 1. In the limit 
that ^ recover our results from the 

main text where scattering is negligible, and decreasing 
the mode volume or increasing the quality factor has the 
same effect on the back-action parameter. In general, 
however, for a given value of ® for an empty resonator, 
there is a maximum achievable rj, 


3 TIME AVERAGED EXPERIENCED 
INTENSITY 

Here we want to derive Eq. (10) of the main text. In 
order to do so, we multiply Eq. (4) with • dt and in¬ 
tegrate both sides over a quarter of an oscillation period: 


^max — 



which occurs at an optimized particle size of 


(S9) 


kr = 


8Q 


(SIO) 


T]{kr) 



FIG. 1: Plot of the back-action parameter as a function of par¬ 
ticle size after considering particle induced scattering losses 
as described by Eq. (S8). Here we take v = 1 and empty- 
resonator quality factors of Q = 10'*, 10®, 10®. 





FIG. 2: Plot of the back-action parameter as a function of par¬ 
ticle size after considering particle induced scattering losses 
as described by Eq. (S8). Here we take v = 1 and empty- 
resonator quality factors of Q = 10®, 10®, 10^. 


dp 


f{x) Wca(w) 


\rix)\ 


2Vr„ 


dt ■ n{x{t)) ■ f{x{t)) 


(Sll) 

We integrate from the classical turning point (where 
the momentum is zero) to the trap minimum (where the 
momentum is maximal) which relates to a quarter of the 
oscillation period T. Using Eq. (9) of the main text, the 
right side of the previous equation is proportional to the 
time averaged experienced intensity and we can formally 
rewrite Eq. Sll as follows: 


{deKp)t — 


4c 


Ta{uj) 


dp 


fix) 

\f'ix)\ 


(S12) 


To proceed we make two approximations: First we ap¬ 
proximate the oscillation period T in the high back-action 
regime as T « 4 ^^* . In particular, the particle moves in 
a square well with length Sx = 2xr = ‘2xt, where Xt is the 
classical turning point, and fmax is the maximum velocity 
in the middle (minimum) of the potential. Additionally, 
in the high back-action regime the particle signihcantly 
changes its momentum only around the classical turn¬ 
ing point Xt when it hits one of the edges of the square 
well. Since the momentum change occurs in a narrow 
region, we can approximate in the integral f{x) « fixt) 
and \ f'ix)\ « \ f'{xt)\. These approximations lead to the 
following equation: 


(Axp) 


t 


cep 2 f{xt) 
a(w) Xt \f'ixt)\ 


(S13) 


where fmax 'Pmax = 2 Ekin, with Akin being the maximal 
kinetic energy in the trap. Now we can normalize this 
time averaged experienced intensity with the time aver¬ 
aged experienced intensity of the optical tweezer regime. 
We begin with expanding the potential Eq. ( 8 ) for small 
W- 

Ut^x) = -2hf^El^^J{x). (S14) 


Since f{x) only varies between 0 and I, it follows that 
the trap depth 5Ut is given by 

5UT = 2h—El —(S15) 
K ° 1 + A2 ^ ^ 


Next we insert Eq. (7) of the main text, and assume 
that the particle is tightly trapped [kSx ^ 1 ) around 
the point of maximum intensity. As the change in intra¬ 
cavity photon number is negligible, we can approximate 
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n(xp) « n{xQ). Eq. (9) from the main text then predicts 
that 

(4xp,T)t = ^^^^^n{xo)f{xo) (S16) 

in the optical tweezer regime. For the case where the par¬ 
ticle is trapped around the antinode of the fundamental 
mode of a Fabry-Perot cavity: 

(/exp.T)i = -^6Ut. (S17) 

Normalizing Eq. S13 with Eq. S17 and ensuring that 
SU = SUt for all 77 yields: 

{Iexp)t ^ 2 fjxt) Ekin 7S1S'l 

(W)* l/'(^*)l 5U ■ 

Surprisingly this equation is valid for all 77 as long as 
the particle is confined sufficiently close to the antinode. 
Fig. S3 shows the excellent agreement between the nu¬ 
merical simulation and the analytic solution obtained by 
Eq. (S18). Taking the limit 77 —)■ 00 of Eq. (S13) implies 

{^exp{v))t 


(texp .T )t 



FIG. 3: Plot of the time-averaged experienced intensity of 
the particle as a function of back-action parameter {/exp(7;))t, 
normalized by the value in the optical tweezer regime 77 ^ 1 . 
The two cases are set to have equal trap depth. The plot is 
numerically calculated for the case of trapping in the funda¬ 
mental mode of a Fabry-Perot cavity f{x) = cos^{kx) with 
resonant positions kxr = fo ■ The red line shows the excellent 
agreement between Eq. (S17) and the numerical simulation of 
Eq. (9) in the main text. 


Xt —>■ Xr which reproduces Eq. (10) of the main text: 


lim (/exp(??))t 

7 ]—>00 


fi^r) .^kin 
Q;(a;) \f'{xr)\ Xr 


(S19) 


4 OPTIMIZATION OF THE HARMONIC 
BACK-ACTION REGIME 


Here we want to maximize the spring constant 
kopt = fchb + kr given by Eq. (11) in the main text. 


^hb = ii^o) describes the first term in 

Eq. (11) and originates from changes of photon number 
with particle position, whereas kx is the familiar term 
known from optical tweezers. The optimization is done 
for a fixed experienced intensity {Iexp)t if we consider 
two trapping modes of a cavity. As a result we will 
derive Eq. (12) and Eq. (13) in the main text and 
conclude how to optimally choose the laser detunings for 
the trapping modes. 


We focus on the regime where the trap minimum 
xq is located roughly at a distance ~ ^ away from both 
resonant positions, where the photon number n{xp) can 
be linearized around the trap minimum xq for each trap¬ 
ping mode i: ni{x) « ni{xQ) + n[{x[))(x — Xq). A linear 
change in photon number with displacement implies a 
harmonic trap, because the force is proportional to the 
photon number (see Eq. (4) and Eq. ( 6 ) in the main text 
and note that f-{x) « /((a^o) ~ for kSx 1 ). 

Using Eq. (9) in the main text, the term proportional to 
7 i'(a:o) does not contribute to the time-averaged intensity 
due to the harmonic motion. Under these circumstances 
Eq. (S16) is valid in the harmonic back-action regime 
as well and the particle experiences the following time 
averaged intensity from each trapping mode i: 

(T .\ 

exp.7 t KWm 1 + (77/((x„))2(a;o - Xri)^’ 

(S20) 

where we linearized the mode profiles around their reso¬ 
nant positions in Eq. (7) in the main text. This is a good 
approximation if the the width of the intensity peaks is 
smaller than the spatial variations of the mode profiles, 
which is the case for 77 ^ 1. Now can write the contri¬ 
butions to the first term fchb of Eq. (11) in the main text 
as: 


fchb.« ° — 2 - 

{.l + {vf[{x.^)Y{xo-x,,Yf 
(S21) 

Expressing the optical tweezer term kx in the same way, 
we can write kopt,i in terms of {Iexp,i)t- 


_a{w) 1 

^opt,2 — \-^exp4/t 


ceo 


fi(xo) [1 +F 


2ri 


- fiixo) 


(S22) 

Xi = \xiifi{xri){x-ci — Xq)] physically describes the ratio 
between half of the width of an intensity peak ^ ^ 

and the distance of the respective resonant position of 
mode i from the trap minimum — xqI- The spring 
constant is maximized for = 1 for which Eq. (S23) 
reduces back to Eq. (12) in the main text. For 77 , >> 1 
and Ti = 1 the contribution to the spring constant 
proportional to f"{ xq) can be neglected and the spring 
constant purely arises from changes of photon numbers 
with particle position. In contrast, for 77 <?; 1 we can 
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neglect the contribution proportional to reaching 

again the optical tweezer regime. 

Eq. (13) in the main text is derived by forming 
the ratio of these two contributions to the spring con¬ 
stant ^ and comparing the two experienced intensities 
necessary to create the same spring constant in each 
regime. To derive this, we also assume that the trapping 
modes consist of the first and second modes of a Fabry- 
Perot cavity, which have equal back-action parameters 
? 7 j. We also use that (/exp,i)t/ 2 (a;o) « {Ie^p, 2 )tfi{xo) 
using Eq. (4) in the main text with Eq. (S16) and 
|/i(a^o)l ~ 1/2 ( 2 ^ 0 ) I close to the trap minimum. 


5 POSITION DEPENDENCE OF SCATTERING 
RATE IN THE HARMONIC BACK-ACTION 
REGIME 


The explicit position dependence of the scattering rate 
«^scat(a^p) contributes positively to the optical spring con¬ 
stant in the harmonic back-action regime with an addi¬ 


tional term proportional to 

(S23) 

where now k = Kex +Kint +Ksca.tixo) and the back-action 
parameter is reduced by particle scattering 77 = ^. 

Thus in the harmonic back-action regime we can de¬ 
fine an effective back-action parameter which is the sum 
of the two terms in the brackets. Interestingly we can 
choose with laser frequency ri = 0 , which sets the term 
coming from frequency shifts to zero. This makes the sec¬ 
ond term arising from position dependent particle scat¬ 
tering experimentally accessible, which would be a factor 
~ (kr)^ smaller than the term arising from frequency 
shifts in its optimized case (ri = 1 ). 


^hb.i — \-^exn,i/1 /. / \ 

ceo Ji(xo) 
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